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(<•-») . Abstract. The return map for planar vector fields with nilpotent 

£SJ ' linear part (having a center or a focus and under an assumption 

generically satisfied) is found as a convergent power series whose 

, ~3 n ' terms can be calculated iteratively. The first nontrivial coefficient 

. is the value of an Abelian integral, and the following ones are 

' explicitly given as iterated integrals built with algebraic functions. 

+-> . 

d . 1. Introduction 

i__i: The study of planar vector fields has been the subject of intense in- 

vestigation, due to their importance in applications, and in connection 
j> ■ to Hilbert's 16th Problem [21]. Significant progress has been made in 

NO : the geometric theory of these fields, as well as in bifurcation theory, 

^ : normal forms, foliations, and the study of Abelian integrals (see the 

CO ' recent book [7]). 

■ First return maps have been studied in relation to existence of closed 

orbits; more generally, return maps are important in a large array 
q : of applications (see [13] and references therein) and also in logic, in 

connection to o-minimality |19j . 

A fundamental result regarding the asymptotic form of return maps 
states that if the singular points of a C°° vector field are algebraically 
isolated, there exists a semitransversal arc such that the return map 
admits an asymptotic expansion is positive powers of x and logs (with 
the first term linear), or has its principal part a finite composition of 
powers and exponentials [13], |20j. 

In the case when the linear part of the vector field has non-zero eigen- 
values there is a good understanding of the return map [2] , [3] , [3] , [S] , [S] , 
[S], as well as for perturbations of Hamiltonians [TU], [TT], or for per- 
turbations of integrable systems [12]. In the general setting, however, 
there are few results available [I], [6], [17], [T8] . 

The present paper establishes an iterative procedure for calculating 
the return map as an integer power series for generic vector fields with 
nilpotent linear part in the case of a center or a focus. The first few 
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coefficients in these series are explicitly given. Using an algorithm given 
here these coefficients can be found explicitly up to any (finite) order 
in terms of iterated integrals involving algebraic functions. 

A fundamental result concerning planar vector fields analytic near a 
stationary point (0, 0) states that if the linear part is not zero, having 
both eigenvalues zero, then such fields have the normal form near (0, 0) 

(1) x = y, y = a(x) + yb(x) 

with a(x),b(x) analytic at 0, with a (respectively, b) having a zero of 
order at least two (respectively, one) at x = [T5] . 

Furthermore, if the origin is a center or a focus, and if the linear part 
of this system is not radial (i.e. does not have the form x = Ax, y = Xy), 
then after an analytic change of variables (CQ) can be written as [16] 

(2) z = -wf(z) + z l+1 g(z), w = kz 2k ~ l f{z)+kwz l g{z) 
where 

(3) 1 < k < I + 1 and /(0) ^ 

The present paper studies the first return map for ([2]) under the sup- 
plementary assumption: 

(4) k^l + 1 

2. Main result 
The main result is the following: 
Theorem 1. 

Consider the system |D) with f(z),g(z) analytic at and satisfying 

(f^]) ; and its transversal T : f(z) w = z l+1 g(z) (z > 0, small). 
Then the first return to T of the solution with z(0) = e (with e small 
enough) is at the point with z = Z(e) analytic in e, with 

Z(e) = e+ enZ " 

n>l-k+2 

where Z n can be calculated iteratively in terms of iterated integrals in- 
volving algebraic functions. 

The proof of Theorem[T] shows that T is indeed transversal. Also, 
the proof contains, and it is built around, an algorithm for calculating 
iteratively the coefficients Z n ; in particular, the first two of them are 
calculated: fUD, ([62} with the notations (|9|h (fT]), (|271), (|29|), (|52|). 

Necessary conditions for a system to have closed orbits (thus a cen- 
ter) are discussed in §3.61 
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3. Proof of TheoremCQ 

3.1. Normalization. Eliminating the time in (j2J) we obtain 

dw kz 2 ^ 1 f(z) + kwz l g(z) 
^ ~dz~ = -wf{z) + z l+l g(z) 

Solutions of (|2J) provide smooth parametrizations for solutions of (jSJ). 
There are points where the graph of a solution w(z) of flSJ) has vertical 
tangents; at such a point the parametrization (J2J) prescribes that w(z) 
is continued with another solution through the same point, and having 
a vertical tangent as well. 

It is convenient to straighten the curve T using a substitution: let 

(6) w = z l+1 F(z) +wi 

with 

(?) 

(note that F(z) is analytic at since /(0) ^ 0). Equation (jSJ) becomes 
(8) 

dz 

with the notations 



(9) p = Z-ib + l, p>l 

(10) = l + ^ 2p F 2 (^) = A F (z) 

dF 

(11) = 2z— + 2{p + 2k)F{z) = B F {z) 

dz 

The problem translates into the study of the first return to the posi- 
tive z-axis of solutions of (jSJ) which start close enough to the origin, i.e. 
which satisfy the initial condition w(e) = with e > small enough. 

It is convenient to introduce the small parameter e in the equation. 
With the notation 

(12) z = ex, wi = e k y 
equation (jSJ) becomes 

(13) = -2kx 2k - x A{tx) - eVyx^- 1 B(ex) 

While the initial condition W\{e) = becomes 7/(1) = 0, it is useful 
to study solutions with the more general initial condition 7/(77) = with 
77 in a neighborhood of 1. 
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It is assumed that 

(14) rj G (1 — r/ , 1 + ?7 ) for some fixed r/o G (0, 1) 

3.2. General behavior of solutions of (1131) . In parametrized form 
equation (TT3"j) is 

(15) x = -2y, y = 2k x 2k ~ l A{ex) + e p yx v+k ~ x B{ex) 

Consider the solution of f|T5|) with the initial condition x(0) = r], y(0) = 
for some 77 > 0. Let (x ,y ) be the solution of the system (TT5T) for 
e = 0: 

x = -2y , y = 2kx 2k ~\ x (0) = 77, y (0) = 
This solution satisfies 

(16) yl + 4 k = v 2k 

and clearly 

(17) x(t)=x Q (t) + 0(e), y(t)=y (t)+O(e) (e - 0) 
for any finite interval of time. 

3.3. Positive solutions of (I13h for x > 0. Lemma[2]shows that there 
exists a unique solution of (TT3"]) so that 2/(77) = and y > for e = 0. 
It also shows that this solution is defined for x G [0, 77] and establishes 
an iterative procedure for calculating its power series expansion in e. 

With the substitution y = u 1 ^ 2 equation (fl~3"j) becomes 

(18) — = -2k x 2k - l A(ex) - eV/ 2 z^" 1 B(ec) 
ax 

Lemma 2. There exists > so that so that the following holds. 

Let rj with [T$ ). For any e with |e| < e equation / TJgj) loirj/i /j/ie 
condition 14(77) = /ias a unique solution u = u(x;e, 77) /or x G [0,77]. 
VKe /iaf e w(a;; 0, 77) > /or x G [0, 77), and u(x; e, 77) zs analytic in e and 
77. 

Proof of LemmalB 
Denote 

(19) P(z; 77, e) = / \-2kt 2k - x A(et)] dt 

7] - x J n 

2k— 1 pi 

''" / 2ka 2k - 1 A (677(7) rfa 

J xln 



Heuristics. A particular vector field was studied in [6], unveiling 
the main ideas involved. In the present general case a few heuristic 
considerations are the following. 
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Looking for solutions of ( {TBI with y{rj) = we obtain that 

y(x) ~c(r]-x) 1/2 (l + o(l)) (x -> rj) with c 2 = 2kr ] 2k - 1 A(er]) 

Note that v4(e^) = 1 + 0(e) > for e small enough, in view of ffTOj) . 
Therefore 

(20) u = y 2 = (r] — x)P(x; rj, e) + A(x) 

with A = o(r] — x) as x — > rj. In fact, substitution of (1201 in ffl8l) gives 

(21) A ~ const (ji - xf /2 

These observations motivate the following substitutions. 
Denote 

(with the usual branch of the square root for x/r] < 1). Let 

(22) 5 = er] 

Note that (TT9T) can be written P(x; 77, e) = r] 2k ~ 1 P(^; S) where 

(23) P(£; 5) = r 2 / 2A; a 2 *" 1 A (5a) da 

Ji-e 



2k [\l -^sf^A^l-^s)) ds 
Jo 



With the substitution 

(24) u (x) = ( V - x)r, 2k - 1 [P(£; 5) + v (£; 5)] 
equation ({18]) becomes 

(25) £ ^ + 2 v = 2P £(1 - e 2 ) p+fc ^ (P + ^) 1/2 5 (5(1 - e 2 )) 
Note that y(r/) = implies = and then necessarily v(0) = 

by (USD, P5, (EU). 

LemmaH follows with eo = <$o/(l — ^o) if we show the following: 

Lemma 3. These exists 5q > so that for any 5 with \S\ < Sq equation 
( fl3]) has a unique solution v = v (£; S) so that v(0) = and this solution 
is defined for £ 6 [0, 1]. 

Moreover, v (£; •) is analytic for 5 G C |o| < <5 . 

T/ie terms of its power series 

(26) „(£;£) =^^ n (£) 
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can be calculated recursively. In particular, the first terms are as fol- 
lows. 

Using the notation^ 



(27) 



(28) * fc (0 = c 



and 



ds 



x _p (1 - S 2fc )V2 



daa k (l-a 2k ) 1 ' 2 



(29) B = 2(p+2k)F(0) = B F . , B x = 2 (p+2k + l)F'(0) = 
we /mve: 

(30) = So^p,fe(0 
and, moreover, for p > 2 we have 

(31) 5) = 5 p v p (0 + 5 P+1 B 1 % +k+1 , k (0 + O (6 P+2 ) 
while for p = 1 we have 

1 



(32) v(Z;5) = 5v 1 (£) + 5 i 



+ O (<5 3 ) 



Proof of Lemma\^ 

Multiplying (1251) by £ and integrating we obtain: 



-2 



C + 25 p t t 2 (l - t 2 )^- 1 [P(f, 5) - v(t; 5)] 1/2 B (5(1 - t 2 )) dt 
Jo 

Since v (0) = then the constant C must vanish. It follows that v is a 
fixed point (v = J[v\) for the operator 

(33) 

J[v] (£;6) = 25 P C 2 [ t 2 {l-t 2 Y +k - 1 [P{t;5)+v{t;5)} 1/2 B(5(l-t 2 )) dt 
Jo 

= 28 p £ [ s 2 (l-£ 2 s 2 ) p+k - 1 [P(£s; 5) + v(£s; 5)} 1/2 B (6(1 - £ 2 s 2 )) ds 
Jo 

Estimates for P(t; 5): 

Let r > be small enough so that F is analytic in a neighborhood 
of \z\ < r, and so that \A(z) — 1\ < cq for all \z\ < r for some Co with 
< c < l/(2k) (see ([TUI)). 

Then from (|23|) . for £ G [0, 1] and |<5| < r we have the upper estimate 

(34) \P(t,S)\ <2k(l + c ) 



1 ^ > p,/c(C) is expressible in terms of the incomplete beta function. 
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and the lower estimate 
(35) 

\P(Z;S)\>2k I {l-esY k - x ds-\2k !\l-es) 2k - 1 [A(6(l-es))-l] ds\ 
Jo Jo 

>2k f (i - s) 2 *- 1 ds - 2k I {l-^sf k ' 1 \A(8{l-i 2 s))-l\ds 
Jo Jo 

> 1 - 2kc > 

The operator J is contractive: 

Let M = sup| 2 |< r \B(z)\. 

Let fi be a number with < /i < 1 — 2kc§. 

Let 5 > be small enough, so that 

(36) «|[2 M 1 + C „) + ,]^M<,, % 3{1 _ 2 ^_^ /2 <1 

Let 23 be the Banach space of functions /(£;<5) continuous for £ G 
[0, 1] and analytic on the (complex) disk \5\ < <5 , continuous on \5\ < 
So, with the norm 

||/||= sup sup \f(Z;5)\ 

£e[0,l] !<5|<<5o 

Let B„ be the ball B^ = {f eB; \\f\\< fi}. 

The operator J defined by ( 133]) is defined on B^. Indeed, for / G B^ 
we have 

(37) \P(t;5) + f(t;5)\ > \P(t;5)\ - \f(t;S)\ > 1 - 2kc - /x > 

1 /2 

therefore, since / and P are analytic in S, then so is (P + /) , and 
therefore so is J[f}- 

Also, J By, C ^ since for f e B^, using (|53p, (|5B|) we have 

|W6<*)| < W^[2k{l + c,)+^ 2 M<^ 

Moreover, the operator J is a contraction on B^. Indeed, using the 
estimate 

\(p+h) 1/2 -(p+f2) l/2 \ < \h-M I sup \p+fr 1/2 < l/j /21 



2/ 6 4' J ' -2(1-2*^-^)1/ 



2 



(by fl37|) ) we obtain 



| JU -Jh\< 4h - h\\ with c = Kw -^—^ < , 



(by m 

Therefore the operator J has a unique fixed point in B^, which is 
the solution v(£;6), analytic for 5 G C with \8\ < 5 . 
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A recursive algorithm for calculating the power series in e: 
To obtain the power series PBjl substitute an expansion v(£;5) = 
+ En>o ^" u n(0 m (j25|) - It follows that for n < p we have £v' n + 2v n = 
with v n (0) = 0, therefore w n (0 = for n < p. 

Using (1231 . ffTUl) we obtain for P(£; 5) a power series in S, with coef- 
ficients polynomials in £: 

(38) P(£; 5) = P (0 + E <5 2p+m ^WO 

m>0 

where 

1 _ (l-$ 2 ) 2k 

(39) P (0 = { e ^ ' 

and 

2A; 1 — (1 — ^ 2k+2 P +m 

(40) pjp+m (0 = P2,m 7. o; 

2p + 2k + m 4 Z 

with the notation 

(41) F\z) = Y. F ^ m 
and in particular 

(42) P 2i0 = P(0) 2 

Substitution of (1261) . (|38|) . followed by expansion in power series in S 
give 

(43) (P(£; 5) + 5)) 1/2 P (5(1 - £ 2 )) = £ 

n>0 

where P n = P„[w p , . . . , v n , £]. 

From (|25p and (j4"3"l) we obtain the recursive system 

with the only solution with f„(0) = given recursively by 

(44) v n (C) = 2C 2 t t\l - i 2 )^- 1 R n _ p (t) dt 

Jo 

The first terms: 

To calculate the first few R n note that 

(45) {P&5)+v&5)) 1 > 2 B(5{l-e)) 

= b p (o i/2 + 5PiP (o l/2 (i - a + 5P JSp; + (* 2 ) 
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where B , B\ are coefficients in the expansion B(z) = B + zBi+O(z 2 ), 
and in view of ffTTj) . they are (1211 . 
In particular 

(46) v p (o = 2B r 2 1 t\i - ey +k - x p Q {t) x i 2 <tt 

Jo 

which equals fl30l) . 

For p > 2 relation (j45]) is 
(47) 

(P + t;) 1 / 2 5 (5(1 - £ 2 )) = ^(O 172 + ^iPo(0 1/2 (l " ^) + O (5 2 ) 

and using fT4"5l) we obtain fl3Tj) . 

For p = 1 we have 0(5 2p ) = 0(5 P+1 ) and there is one more term in 
the second nontrivial coefficient of i>(£; 5). The calculation is straight- 
forward: using relation ( )45i) for p = 1 and (jSJ) we obtain ( 1321) . 

□ 

The following Corollary gathers the conclusions of the present sec- 
tion. Many quantities depend on the function F, see ([7]), (flOl) . (jTT|) . 
(|23|) . and we add the subscript F for them: 

Corollary 4. There exists €q > so that for any rj with (Lfy and e 
with |e| < e equation l[TS\) has a unique solution y(x) on [0, rf\ satisfying 
y(r)) = and y > on [0, rj) for e = 0. 

Moreover, this solution has the form 
(48) 

y{x) = <f) F (x; e, rj) = (r]~x) 1/2 f] k - 1/2 P F ^1 - ^; er^ + w F ^1 - ^; e 

rat/i -Pf(£; 5) = 5) gwen 6?/ ( Tggj) andv F ((; 8) solution of (25)) with 
v F {0;5) = andv F {£;0) = 0. 

The map (e, 77) 1— > t>F(£; e^) analytic for \e\ < €q and 77 as in ( [I^P - 

□ 

3.4. Solutions of (1131) in other quadrants and matching. 

3.4.1. Solutions in the four quadrants. CorollaryH] gives an expression 
for the solution y(x) of (fT3|) for x > and y > 0. Solutions in the other 
quadrants are found in the following way. 

Let e with |e| < e with eo given by LemmaO 
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(i) Let yi(x) = 0^(x;e, 77) the solution ( 1481) of ( TT3l) . defined for x G 
[0, 77], with t/i(??) = 0. We have 7/1(2) = (r] 2k — x 2k ) l l 2 + 0(e) in view of 
f fl6|) . ffT7|) . and we will refer to 7/1 as a "solution in the first quadrant". 

Solutions "in the other quadrants" are obtained as follows. 
For a function F(z) denote by JiF the following functions: 

(49) (J 2 F)(x) = (-iy+ k F(-x), (J 3 F)(x,y) = (-lT +k - 1 F(-x), 

(J 4 F)(x) = -F(x) 

Note that J2J3 = J4 and Aj 2 p = Ap. 

(ii) It is easy to check that the function 7/2(2) = 0j 2 f(— 2; e, 77) is 
a solution of ( TT3l) . It is obviously defined for 2 G [—77,0]; we have 
y 2 (-rj) = and y 2 (x) = [rf k - x 2k ) 1/2 + 0(e). 

(iii) Similarly, the function yz{%) = — 0j 3 f( - 2; v) is a solution 
of f[TH]) defined for x G [—77,0]. We have 7/3 (— 77) = and 7/3(2) = 

(iv) The function 7/4(2;) = — </>j 4 f(2; e, 77) (a;; e, 77) is a solution of ([TBI) 
defined for x G [0, 77]; we have 7/4(77) = and 7/4(2) = —(rf k — a; 2 * 1 ) 1 / 2 + 
O(e). 

3.4.2. Matching at the positive y-axis. Let 77, 77 satisfying (fT4"|) and let 
7/1(2) = (Pf(x; e,rj) solution as in (i), for x G [0,7/] and 7/2(2) = 
(Pj 2 f(—x; e, 77) solution as in (ii), for a; G [—77, 0]. 

The following Lemma finds fj so that 7/1 (0) = 7/2(0), therefore so that 
7/1 is the continuation of t/ 2 : 

Lemma 5. Let |e| < e . £e£ 77 so that \rj — 1| < Ci|e| 7x>rf/i Ci small 
enough so that Cie < 7/ . 

There exists a unique fj = 7/ + O(e) so t/iai 

(50) F (O;e,7/) = j2F (O;e,7/) 
Denote this 

(51) ^ = A/>(7/,e) 

Moreover, A/ir(7/,e) depends analytically on 77 and e /or |e| < ei /or 
ei small enough. We have \fj — 1| < C2 1 e| /or some c 2 > 0. 

Furthermore we have, in the notations (27\ ), ( [HP , §E3\) , and with 



(52) 



o P = (-1) 
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that for p > 2 



(53) fl = M F ^e)=r 1 + eV^ l ^^{l + e p ) 



and /or p 



(54) r/ = Ar F fa6)=r / + 6r / 2 ^ ^ (1) (1 + g t ) 

2 a O^gQ^jMl) + (l + gi)/fcB 2 $ 1|fc (l) 2 3 ^ 
77 2fc 

Proof of Lemma{5[ 

Using (|48p equation (1501) is equivalent to solving the implicit equation 

(55) Gifi, rj, e) = where 

G(fj, V , e) = fj 2k [P JaF (l; efj) + v. hF (l; efj)] - r, 2k [P F {\- eq) + v F (l; e V )] 

which is a function analytic in (fj, r], e) by Lemma[21 
We have G(r],n,0) = (by (EH]), (ESJ), (EH) and 

therefore by the implicit function theorem and using the compactness of 
the interval \rj— 1| < cieo, equation G(fj, 77, e) = determines 77 = 77(77, e) 
as an analytic function if |e| < E\ for t\ small enough. 
Since 77(77, 0) = r\ we have, for |e| < E\ and (77 — 1| < c\e\, 

I77C77, e) — 77I < |e| sup — =ci|e| 

\e\<ei t \v-l\<cie! Ot 

therefore 

\fj(v>e) ~ !| < l r ?( 7 ?7 e ) - ^1 + 1*7 ~ 1| < ( c i + c 'i)l e l = c 2|e| 

We can assume c 2 ei < 770 by lowering e%. 

The expansion of 77(77, e) in power series of e is found by introducing 
the expansions ([38]), ([26]), (]4l| in ([55]) and noting that P (l) = 1, 
P J2F;2 (1) = p 2 (i) (see dSSJ), dH), and that the coefficients of 

5j 2F in m are /3 J2F;0 = (-l) p+fc A) and Bj 2F;1 = (-l)^ 1 ^. 

□ 
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3.4.3. Matching at the negative y-axis. Let rj, fj satisfying (115H and 
1 77 — 1| < ci|e| as in Lemma[5j Let fj given by LemmaO Consider 
the solution 2/3(2) = — 4>j 3 f(— x\ e, 77) as in (iii), for x G [—77,0], with 
y~z{— fj) = 0- Therefore 7/3 is the continuation of 7/2- 

Let 7/4(2) = — 0j 4 f(x; e, 77) be a solution of (fl~3l) as in (iv), for a; G 

The following Lemma finds r\ so that 7/3(0) = 7/4(0), therefore so that 
7/4 is the continuation of 7/3: 

Lemma 6. Lei |e| < e\ and \fj — 1| < C2 1 e | 7tra£/i ei small enough so that 
fj satisfies (T73]). 

There exists a unique fj > so that 

(56) - 0j 4 f(O; e, fj) = -0j 3 f(O; e, 77) 
Moreover, we have 

(57) f}=Mj 3F {fj,e) 

where Nf denotes the function (5l\) of Lemma{5{ 

Therefore fj depends analytically on fj and e for \e\ < £2 for £2 small 
enough. We have \fj — 1| < C3 1 e| for some C3 > 0. 

Furthermore, the first coefficients of the e series of J\fj 3 F{rj,e) and 
Nf{j], e) coincide: 

(58) Ar j3 F(v,e)=Ar F (v,e) + 0(e^ 2 ) 
Proof. 

Note that the equation (I56I) for fj = 7/(77, e) is the same as the equa- 
tion (|55|) for fj = 77(77, e), only with F replaced by J 3 F. Noting that 
^j 3 F;2(l) = ^2(1) and that Bj 3F;0 = 9 P B , B JzF] i = ~0 P B 1 Lemmal 
follows from Lemma[5j □ 

3.5. The first return map. Let 77 satisfying (fl4|) and 69 as in Lemmal6l 
Then 77 given by Lemma[6] is the first return to the positive of 
the solution with x(0) = 77, 7/(0) = and it is analytic in e and 77: 

(59) 7] = A/" j3 F(A/>(77, e ), e ) 

At this point it is enough to take 77 = 1. Going back through the sub- 
stitutions (fl2l) . ([6]), we obtain that the first return map to the transver- 
sal w = z l+1 F(z) (with z > small) of the solution of ([2]) with z(0) = e, 
w(0) = e l+1 F(e) is attained for 



(60) z = eMj 3F {M F {l,e),e) = e + ^ Z n e 

n>p+l 



FIRST RETURN MAP 



13 



and the terms of this convergent power series in e can be calculated 
recursively. 

In particular, the first terms are obtained from (1351) . (154|) . (|58|) : for 
p > 2 we have 



(61) z = e + 2e P +l g °y ) (l + fl p ) 



+ 



2A: 



and for p — 1 



(62) z = e + 2e 



+ 2 e J 



2 go^fcO) 

2A; 



^-gi^grg^O) + (l + gi)/fcBg^ 1|fc (l) ; 
2A; 



+ 0(6 3 ) 



3.6. Closed trajectories. Relation (1601) together with the construc- 
tive method presented allows to calculate, recursively, the coefficients 
Z n , and therefore to decide whether the origin is a center or a focus: 
the origin is a center if and only if all Z n vanish. 

As a practical matter, for each value of p = I — k + 1 all the power 
series should first be properly ordered (the order in which the terms 
appear does depend on the value of p) and then the series can be 
calculated (at least in principle) to any order no- The conditions Z n = 
for n < no can be written in terms of the function F, and they are 
necessary conditions for the fixed point to be a center. 

The first such conditions follow from the first two nontrivial terms 
calculated here: using (IBTl) . (16"2"1) for the origin to be a center we must 



have Z 



P+i 



J p+2 



which implies B (1 + 9 P ) = and B 1 (1 - 9 P ) 



0, which in turn means that if p + k — 1 is odd then we must have 
F'(0) = 0, while if p + k — 1 is even, then we must have F(0) = 0. 
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